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Abstract 

The dynamics of the one-dimensional random transverse Ising model with both nearest-neighbor 
(NN) and next-nearest-neighbor (NNN) interactions is studied in the high-temperature limit by 
the method of recurrence relations. Both the time-dependent transverse correlation function and 
the corresponding spectral density are calculated for two typical disordered states. We find that 
for the bimodal disorder the dynamics of the system undergoes a crossover from a collective-mode 
behavior to a central-peak one and for the Gaussian disorder the dynamics is complex. For both 
cases, it is found that the central-peak behavior becomes more obvious and the collective-mode 
behavior becomes weaker as Ki increase, especially when Ki > Jj/2 (J, and Ki are exchange 
couplings of the NN and NNN interactions, respectively). However, the effects are small when the 
NNN interactions are weak (Ki < Ji/2). 
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I. INTRODUCTION 



There has been a considerable interest in the study of the dynamics of quantum spin 



systems in the past few decades 



2|, |3j, and the calculation of dynamic correlation functions 



remains a highly nontrivial and real challenging task until now. Some exact results have been 
obtained for the one-dimensional (1-D) pure quantum spin models, e.g., the 1-D transverse 
Ising model and the 1-D XY modelj^, 5, (|. Recently, more attention has been paid to 
the investigation of the dynamical behavior of disordered systems^, 



11 



12 1, which 



can be applied to describe the properties of many materials such as window glass, magnets 
with frozen-in disorder, etc.. One simple but important example of such systems is the 1-D 
random transverse Ising model (RTIM). 

The dynamical behavior of the 1-D RTIM with the bimodal distribution is studied by 
Florencio and Barrotojsj], and it is found that the dynamics undergoes a crossover from a 
central peak behavior onto a collective mode one. Recently, we have investigated the effects 
of Gaussian disorder on the dynamics of the 1-D RTIM 12], and have found that there are 
two crossovers when the standard deviation of random variables is small and there is no 
crossover if the value of the standard deviation is large enough. Besides, the dynamical 
behavior of the random-bond transverse Ising model with four-spin interactions 13j and the 



disordered XY chain 



B y 



ll| have been studied. 



In the work mentioned above only nearest-neighbor (NN) interactions are considered. 
To our knowledge, no related results of disordered quantum spin systems with next-nearest- 
neighbor (NNN) interactions have been reported so far. However, the work of Sen has shown 
the role of second neighbor interactions on the relaxation in s = 1/2 pure quantum spin 
chains 1J]. The results show that the dynamical correlation in the NNN transverse Ising 
chain is noticeably different with that of the exactly solvable NN transverse Ising chain. 
Therefore, it is expected that considering additional spin-spin interactions in disordered 
systems will make some differences in the dynamical process. 

Our main interest is to investigate the effects of NNN interactions on the dynamics of 
the 1-D disordered quantum spin systems. It is well known that the interactions between 
spins may be complex in real materials. For studying the properties of real systems theoret- 
ically and experimentally, the easy way is to consider a model with only the dominant NN 
interactions. In this paper, we consider the 1-D RTIM with both NN and NNN interactions. 
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The 1-D RTIM can be used to describe the dynamical property of many condensed-matter 



systems like the quasi-one-dimensional ferroelectric crystals (e.g., PbHi_ x D :r Po4)ll9|, ]20|, 
and the Ising spin glass LiHo .i67Yo.833F4 21]. We investigate the cases that the exchange 
couplings between spins or transverse fields independently satisfy the bimodal distribution 
and the Gaussian distribution, respectively. Our calculations are based on the method 
of recurrence relations 



22 



turn many-body dynamics 



231 which is very powerful in the study of classical and quan- 



171 ] . Meanwhile, we also used some reliable approximation 



25] and the Pade approximants. 



schemes such as the so-called Gaussian terminator |3|, ]24], 
It is found in both disorder that the central-peak behavior becomes more obvious and the 
collective-mode behavior becomes weaker when Ki > Jj/2. We also find that the dynamics 
of the system is not sensitive to the property of the NNN interaction whether it is ferromag- 
netic or antiferromagnetic. 

This paper is arranged as follows. In Sec. [IT] we introduce the model used in this paper 
and the method of recurrence relations. Sees. II I II and IIVI give the dynamical results for the 
bimodal disorder and the Gaussian disorder, respectively. Sec. IVl provides conclusions. 



II. MODEL AND METHOD 



The Hamiltonian of the 1-D RTIM with both NN and NNN interactions can be written 



as 



where of 



N N 



(i) 



a 



x, y, z) are Pauli matrices at site i, denote the external fields, while Jj and 



Ki are exchange couplings between NN spins and NNN spins, respectively. The periodic 
boundary conditions af +N = erf are assumed in next calculation, where N is the number 
of spins. For simplicity, we assume that Ki = aJi (0 ^ a < 1) and consider Jj and Bi are 
uncorrelated random variables which satisfy the probability distributions p{J%) and p(-Bj), 
respectively. It is obvious that, in the limit Ki — > in Eq. (CQ), this model can be reduced 
to the 1-D RTIM^Q. 

The dynamical behavior of classical or quantum many-body systems is conveniently ex- 
pressed in terms of dynamic correlation functions. In this paper, we are interested in the 
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time-dependent transverse correlation function defined by 



(2) 



where (• • • ) denotes an ensemble average followed by an average over the disorder variables. 
The spectral density $ (u) (u> is the frequency) which is able to be determined directly from 
experiments is defined as the Fourier transformation of the correlation function, 



iuit 



C (t) dt. 



(3) 



The method of recurrence relations has already bee n ap plied to solve a variety of 



body syste ms, such as the classical harmonic chain 13] , the electron gas 22|, |23| , spin 



many- 



systems 15|, lla, |24j and ergodic theory [18| etc., successfully. In the following we will sum- 



marize this method. 

Consider a many-body system defined by a Hamiltonian H. The time evolution of a 
dynamical operator A is described by the Liouville (or Heisenberg) equation of motion 

dA (t) 



dt 



iLA (t) 



where L is the Liouville operator, LA = [H, A] = HA — AH. The solution of Eq. 
be given as the form of the orthogonal expansion 22|] 



(4) 



can 



(5) 



v=0 



where {f v } are an orthogonal set of basis vectors spanning a Hilbert space S, the coefficients 
a v (t) are time dependent functions representing the projection of A(t) onto f v at t. 

In the high-temperature limit T = oo, the inner product which includes both the statis- 
tical and random averages in our system is described as 



I,19|,H0| 



(x,y) = <xYt), 



(6) 



where X and Y are basis vectors defined in S. 

Set fo = A (0), which gives ao (0) = 1 and a u (0) = for v > by Eq. (JHJ). The basis 
vectors f u satisfy the recurrence relation (RRI) 



f v+1 = iLfv + A„/„_i, v > 0, 



(7) 
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where the coefficients, also known as recurrants, are defined as 

A,= ( } fu,f ; ) s iy>\) (8) 

\Ju-l, Jv-l) 

with /_i = and A = 1. Meanwhile, the coefficients a u (t) satisfy a second recurrence 
relation (RRII) 

A u+1 a u+1 (t) = -^P- + du-i (t), v> 0, (9) 

where a_i (t) = 0, and Oo (i) is the time-dependent correlation function (A(t)A(O)). Ob- 
viously, by choosing / = crj the average spin correlation function is just given by 
C if) = {&j (t) <jj (0)) (Eq. (j2J)), which can be written as the form of moment expansion 

k 

2k 



( — 1 



fc=0 

with 

^2* = ^fW, [H,---[H,a?]---]], (10) 

where fj,2k is the 2fcth moment of C (t). The partition function Z =Trl = 2 W equals the 
number of quantum states of the system. Using the first 2u moments, we can calculate the 
correlation function C (t) by constructing Pade approximants. 

By taking the Laplace transformation of the recurrence relation (RRII), one obtains 

A v+1 a v+ i(z) - S V)0 = -za v (z) +a„-i(z) , ^ = 0, 1, 2 - • • , (11) 

where z = e + iu (e > 0) is a variable of the complex plane. Then one can get the continued- 
fraction form 

1 

a (z) = . (12) 

z + 

A 2 

z + 



z + --- 

Furthermore, it is proved that the spectral density $ (u) (Eq. (j3j)) is able to be determined 
directly by Eq. (jl2]h 

$ (w) = limReao (z) . (13) 

Note that A v are the key quantities for calculating the dynamic correlation functions. 

Generally, only a finite number of continued-fraction coefficients can be determined. So 
it is necessary to use a scheme to terminate the continued fraction. The one that serves our 



251 ]. Suppose the first M recurrants 



model best is the so-called Gaussian terminator 
are determined, in this approximation, the others are assumed to be of the form A t 
z/(A M /M), for v> M. 



III. DYNAMICS FOR BIMODAL DISORDER 



After a lengthy calculation, the first eight basis vectors are exactly obtained by Eq. (j7j). 
In the following, we just give the first two of them: 

h = Brf, 

h = (Ax - B]) aj + P>,K, 2 crJ .,o j + !>>,.!, ,rrj v a) • 
BjJtf+tf + BjKtf+tf. (14) 

The squared norms of the basis vectors are given by Eq. (J6j) as follows: 

(/o,/o) = l, 

(/i,/0 = ^, 

(/ 2 ,/ 2 ) = aT-2at^+^+5^+5^+5|^;+5|^. (is) 

Using the above results, we have calculated the first eight coefficients Ai, A 2 , ■ ■ • , and A 8 
exactly, and the A 9 through the assumption /S. v = u(A M /M) approximately. Meanwhile, 
the first 18 moments are obtained and the correlation function C (t) can be determined by 
constructing the Pade approximants. 

Notice that the coefficients (see Eq. (J5])) are even functions of Ki and the correlation 
functions are determined uniquely by the recurrants. Thus, the dynamical property of the 
system is independent of that the NNN interactions are ferromagnetic or antiferromagnetic. 
Actually, the system is in its paramagnetic phase in the high-temperature limit. Next, 
we only consider the case of ferromagnetic NNN interactions (Ki > 0). We calculate two 
typical cases that the random variables satisfy the bimodal distribution and the Gaussian 
distribution, respectively. 

In the following, we assume that the exchange couplings Ji or the transverse fields Bi 
satisfy the bimodal distribution 

N 

P ({Pi}) = TUP* (A - Pa) + (1-P)S (Pi - AO], (16) 
i 
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FIG. 1: Time-dependent correlation functions C (t) and corresponding spectral densities <J? (u) for 
the case of random bond which satisfy the bimodal distribution in which J a = 1.0 and = 0.4. (a) 
and (b) plot of the pure cases that p = and 1. The results in (c) and (d) are the disordered cases 
for p = 0.25 and 0.75. The central-peak behavior becomes more obvious and the collective-mode 
behavior becomes weaker as Ki increase. The black solid line in (a) is monotonic for t < 3. 



where /3, = Jj or B^ p is the concentration of coupling J a or magnetic field B a and takes 
values from to 1. 

We first consider the random band and uniform field model. In this case, without loss 
of generality we set B^ = B — 1, and choose J a = 1.0 and Jf, = 0.4, which have been used 
in Ref. In this assumption, the exchange couplings Jj change from < B (p < 1) to 
Jj = B (p — 1). The transverse correlation functions C (t) and the corresponding spectral 
densities $ (u) are given in Fig. 1 for several values of bond concentration p. In order to 
show better the effects of the NNN interactions on the dynamics of the system, we have 
considered the cases that Ki = 0, Jj/4, Jj/2 and 3Jj/4, respectively. Obviously, for the 
Ki = case, the results are just of the 1-D RTIM studied by Florencio and BarretoQ]. 



7 



From different cases of Ki shown in Fig. 1, we can see that the system shows a collective- 
mode behavior for small values of p (i.e. p = or 0.25), and exhibits a central-peak behavior 
when p = 1. In general, when the external field is small, the spin-spin interactions play an 
important role, thus the central-peak behavior dominates the dynamics of the system, while 
for large B the dynamical behavior is collective-mode one, which is due to the precession 
of spins in the transverse field. This means that the dynamics of the present model with 
bimodal distributions is similar to that of Ref. s|. 

By comparing the curves for the cases that Ki = Ji/4, J{/2, and 3Jj/4 with those of the 
Ki = case (see Fig. 1), we can see that the dynamics has no evident change if the NNN 
interactions are weak (e.g. Ki = Ji/4). However, there are some obvious differences when 
Ki > Jj/2. The dot-dashed curve for the pure case p — 1 when Ki = in Fig. 1(a) describes 
the dynamics for the exactly solvable limit in which Jj = Bi — 1, now the C (t) is a Gaussian 
function^. Meanwhile, the other curves for the cases that p = 1 when Ki = Ji/4, Jj/2, and 
3Jj/4 in Fig. 1(a), respectively, all behave monotonically but exhibit slower decay than for 
the Ki = case, and are not a Gaussian. On the other hand, the lines shown in Fig. 1(a) 
for the pure case p = indicate that the collective-mode behavior becomes weaker as Ki 
increase. The same results can be also obtained from the corresponding spectral densities. 
As shown in Fig. 1(b) the central peak increases, meanwhile, the collective-mode peak 
becomes lower and the width of the spectral line broaden as Ki increase. For the disordered 
case that p = 0.25 (see Figs. 1(c) and (d)), we also find weaker collective-mode behavior 
if the NNN interactions become stronger. However, for the disordered case that p = 0.75 
when Ki = 3Jj/4 (see the black solid line in Fig. 1(c) ), the C (t) decays monotonically to 
zero, and the dynamics of the system is a central-peak behavior which is not as the case 
that p = 0.75 when K j = (the lines for p = 0.75 when Ki = or Jj/4 in Figs. 1(c) and (d) 
show a disordered behavior which is something between the collective-mode behavior and 
the central-peak one). All the above results indicate that the interactions between spins are 
stronger in our system, and that the effects of the NNN interactions on the dynamics of the 
system cannot be neglected. 

We now consider the random field and uniform band model, in which the transverse 
fields satisfy the bimodal distribution p (Bi) and can take the values B a = 0.6 (p = 1) and 
B b — 1.4 (p — 0), while the exchange couplings are constants (Jj = J — 1, Ki — 0, 1/4, 
1/2 or 3/4). This allows the external fields changing from Bi > J to Bi < J as p increases. 
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FIG. 2: Correlation functions and corresponding spectral densities for the random field model, in 
which B a = 0.6 and = 1.4. The system undergoes a crossover from a collective-mode behavior 
to a central-peak one as q increases. The central-peak behavior becomes more obvious and the 
collective- mode behavior becomes weaker as Ki change from to 3/4. 

The results of C (t) and $ (a;) for different values of p are shown in Fig. 2. The curves for 
p = are the pure cases dominated by the stronger field energy. In this case, the system 
is at the collective mode regime. When p = 1 and 0.75, the correlation functions decay 
monotonically, and thus the dynamics is dominated by the central-peak behavior. However, 
for the disordered case p = 0.25, the dynamics of the system is neither central-peak nor 
collective-mode type, but something between them. Hence, for this model, the system also 
undergoes a crossover from a collective-mode behavior to a central-peak one as p increases 
from to 1. 

The same as the above random bond model, it is found that the central-peak behavior 
becomes more obvious and the collective-mode behavior becomes weaker as Ki increase, 
especially when Ki > Ji/2. From the curves for p — 1 in Fig. 2(a) and p = 0.75 in Fig. 2(c), 
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we can find that the C (t) decays more slowly as Ki changing from to 3/4. Meanwhile, 
the magnitude for $ (a>) at u — increases as Ki ^ (see Figs. 2(b) and (d)). The results 
of p = in Fig. 2(a) and (b) indicate that the collective-mode behavior becomes weaker if 
the NNN interactions are stronger, since the oscillatory curves are less damped. 



IV. DYNAMICS FOR GAUSSIAN DISORDER 



In the following, we assume that the exchange couplings Jj or the transverse fields Bi are 



uncorrelated random variables which satisfy the Gaussian distribution 



12 



26] 



p({fl}) = JJ-J^exp 



(A -I) 2 



(17) 



where f3 denotes the mean value of the random variables and 073 is the standard deviation. 
Next, we discuss two different cases that the random-bond and the random- field model, 
respectively. We find that the effects of Ki on the dynamics are not obvious when Ki < 
Ji/2, which is similar as the above case of the bimodal disorder. In the following, we only 
give the results for Ki = Ji/2 and 3Jj/4 in the random-bond model, and for Ki = 3Jj/4 in 
the random-field model. 

For the random-bond model, the exchange couplings J( satisfy the Gaussian distribution 
while the transverse fields Bi are constants. We keep Bi = B = 1 which sets the energy 
scale, and consider that the mean value J varies from to 2 and the standard deviation aj 
changes from 0.3 to 3.0. Both the correlation functions C (t) and the corresponding spectral 
densities $ (u) are shown in Figs. 3 and 4, respectively. The insets to Fig. 4 present the 
first nine recurrants. 

From Figs. 3 and 4 we can see that for the case of small values of oj (e.g., 0.3), there 
are two typical dynamics: the collective-mode behavior and the central-peak behavior. It is 
obvious that the black solid curve for J = in Fig. 3(a) is a damped cosine function, which 
is due to the precession of spins in an external transverse magnetic field 12j. As J increases, 
the system first shows a weak collective-mode behavior for the case of J < B (e.g., J = 0.5), 
then exhibits a central-peak behavior when J > B (e.g., J = 1.5). 

For the case of large <j,j (e.g., 3.0), it is found that the system only shows a central- 
peak behavior, and there is no crossover. In this case, the strong exchange couplings play 
an important role in the dynamics of the system. That is, the spin-spin interactions are 
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FIG. 3: Correlation functions C (t) for the case that the bonds satisfy the Gaussian distribution 
while the external fields are constants. Two typical cases that aj = 0.3 and 3.0 are considered. 
The mean value J varies from to 2. 

dominant in the competition between the spin-spin interactions and the external fields. It 
is also found that, further increasing the NNN interactions will make the curves of C (t) 
for the central-peak behavior decay more slower and the magnitude for $ (u) at u = 
become larger. This all indicate that the dynamical behavior of the system is sensitive to 
the inclusion of the NNN interactions. 

We next discuss the results of the random-field model, in which the transverse fields Bi 
satisfy the Gaussian distribution, while the exchange couplings remain unaltered ( Jj = J = 
1, Ki = 3/4). Let the mean value B varies from to 2. From Fig. 5 we can see that for 
ob = 0.3 the system undergoes a crossover from a central-peak behavior to a collective-mode 
one as B increases. When ob is large enough (ffg = 3.0), the system only shows one type of 



dynamics and there is no crossover. That is a most-disordered state 12j, which is something 



between the central-peak behavior and the collective-mode behavior. By comparing the 
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FIG. 4: Corresponding spectral densities (to) for the same parameters as in Fig. 3. For aj = 3.0 
the central-peak increases as Ki change from Jj/2 to 3Jj/4. For aj = 0.3 the collective-mode 
behavior becomes weaker as Ki increase. 



results with those of Ref. [l2j, we can find that the oscillatory behavior becomes weaker as 
Ki increase. However, different from the effects of the NNN interactions on the dynamics of 
the random-bond model, increasing Ki in this case will make no difference when erg is large 
enough. That is, the dynamics of the system now is dominated by the disordered external 
field. 



V. CONCLUSIONS 

In this paper, we have studied the effects of the NNN interactions on the dynamics for 
the 1-D RTIM in the high-temperature limit. We have considered the cases that the random 
variables satisfy the bimodal distribution and the Gaussian distribution, respectively. It is 
found in both cases that the dynamical property of the present model is similar to that of 
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FIG. 5: Correlation functions and corresponding spectral densities for the random field model when 
Ki = 3Jj/4. The insets to (b) and (d) present the first nine recurrants. Let the mean value B varies 
from to 2. For ob = 0.3 there is a crossover from a central-peak behavior to a collective-mode 
one as B increases. For <jb = 3.0 the system exhibits a most disordered behavior. 

the 1-D RTIM when Ki < Jj/2. However, the central-peak behavior becomes more obvious 
and the collective-mode behavior becomes weaker as the NNN interactions increase (i.e. 
Ki > Jj/2). It is expected that we can get similar results in other disordered quantum spin 
systems, e.g., the XY and the XYZ models. 
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